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Abstract 



On the basis of the two-site polaron problem, which we solve by exact di- 
agonalization, we analyse the spectral properties of polaronic systems in view 
of discerning localized from itinerant polarons and bound polaron pairs from 
an ensemble of single polarons. The corresponding experimental techniques 
for that concern photoemission and inverse photoemission spectroscopy. The 
evolution of the density of states as a function of concentration of charge carri- 
ers and strength of the electron-phonon interaction clearly shows the opening 
up of a gap between single polaronic and bi-polaronic states, in analogy to 
the Hubbard problem for strongly correlated electron systems. In studying 
the details of the intricately linked dynamics of the charge carriers and of the 
molecular deformations which surround them, we find that in general the dy- 
namical delocalization of the charge carriers helps to strengthen the phase co- 
herence for itinerant polaronic states, except for the crossover regime between 
adiabatic and anti-adiabatic small polarons. The crossover between these two 
regimes is triggered by two characteristic time scales: the renormalized elec- 
tron hopping rate and the renormalized vibrational frequency becoming equal. 
This crossover regime is then characterized by temporarily alternating self- 



localization and delocalization of the charge carriers which is accompanied 
by phase slips in the charge and molecular deformation oscillations and ul- 
timately leads to a dephasing between these two dynamical components of 
the polaron problem. We visualize these features by a study of the temporal 
evolution of the charge redistribution and the change in molecular deforma- 
tions. The spectral and dynamical properties of polarons discussed here are 
beyond the applicability of the standard Lang Firsov approach to the polaron 
problem. 
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I. INTRODUCTION 



Renewed interest in the physics of small polarons in the last few years has been largely 
stimulated by the discovery of new materials with exceptional properties such as the High 
Tc cuprates, the nickelates and the manganites showing giant magneto resistance. It is 
believed that small polarons play an essential role in these materials. Recently developed 
experimental techniques such as femto- second optical spectroscopy, EXAFS and pulsed neu- 
tron scattering, ionchanneling experiments and high resolution angle resolved photoemission 
experiments have been employed to study the properties of polaronic systems. Inspite of 
considerable theoretical efforts in the last few years, the underlying physics of the polaron 
problem has remained largely unresolved. 

As concerns the problem of a single polaron in an empty lattice with local coupling of the 
charge to a set of individual non-interacting local lattice deformations (the Holstein model) 
it is known that self-trapping of a charge carrier occurs when 

i) the gain in localization energy Sp outweighs the gain in kinetic energy D/2, i.e., g = 
2ep/D > 1, D denoting the bare electron bandwidth and 

ii) the relative deformation of the lattice which surrounds the charge carrier essentially 
remains confined to the immediate vicinity of the charge carrier i.e., = £p/uJo > 1, ujq 
denoting the bare local phonon frequency. 

A crossover between essentially delocalized quasi-free electrons and self-trapped electrons 
is known from Quantum Monte Carlo simulations to occur in a regime of parameters 
characterized by the two conditions i) and ii) for small polaron formation, i.e. for ~ 1 and 
a ~ 1. This crossover regime in a narrows substantially upon going from the anti-adiabatic 
(7 = t/uo = Aa^ /{zg) <^ 1) to the adiabatic (7 ^ 1) limit with z denoting the coordination 
number. 

The standard theory of small polarons is based on the so called Lang Firsov (LF) trans- 
formation usually followed by approximations treating the fully localized polaron state as 
the starting point of the theory and then introducing the hopping of the electrons by per- 
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turbative means p[. There is a widespread (however erroneous) behef that such a scheme 
is vahd in the extreme strong couphng anti-adiabatic hmit i.e., 7^1 and a ^ 1. It is not 
surprising that such theoretical treatment is even less able to describe polarons close to the 
crossover regime, which represents the physically interesting and realistic situations. 

As concerns the problem of the many polaron system the principal difficulty is to take 
into account the overlap of the lattice deformations surrounding the charge carriers when the 
density of polarons becomes important. It is expected that, due to destructive interference 
of such local deformations, the whole concept of polarons breaks down, the lattice changes 
its structure and as a result the effective electron lattice coupling becomes strongly reduced. 

It is the purpose of this present work to obtain some insight into the physics of small 
polarons for realistic values of the relevant parameters 7 and a and in particular close to 
the crossover regime. Since all known materials containing small polarons consist of highly 
polarizable small clusters {Ti^~^ — Ti^^ pairs in Ti^O-j, V'^'^ — V'^'^ pairs in NaxV205), 
sometimes embedded in a metallic background (such as the 0(4)^~ — Cu{l)^ — 0(4)^~ 
dumbbells in cuprate High superconductors containing chains separating the metalic 
Cu02 layers), it is not only instructive from a pedagogical point of view to study small 
polaron features on the basis of such small units, but such studies may also apply to true 
physically relevant situations. On the basis of a two-site polaron system (a system consisting 
of two adjacent molecules between which the electrons can hop) we shall illustrate the highly 
non-linear physics going on in such a problem. For a ^ 1 and 7^1 (the so called strong 
coupling anti-adiabatic regime or extreme polaron limit) it will be shown that the dynamics 
of the fluctuations of the lattice deformations is driven by the charge fluctuations, while for 
a -C 1 and 7^1 (the so called weak coupling adiabatic regime) the inverse happens. The 
crossover regime is characterized by alternating localization and delocalization of the charge 
carriers as a function of time. In general the amplitude and the frequency of the intrinsic 
lattice vibrations as well as of the intrinsic electron hopping rate are strongly renormalized. 
Such effects may turn out to be vital as concerns a mechanism for the damping of polaronic 
charge carriers; a feature which is inaccessible within the usual LF approach. 
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The paper is organized in the following way: In section 2 we shall present the model and 
its basic physics as well as our method of an exact diagonalization study. In section 3 we 
discuss the dependence of the kinetic energy of the electron on the adiabaticity parameter 
7 and the coupling strength a and examine the limitations of the LF approach. Section 4 
will be devoted to the study of the one particle spectral function in view of distinguishing 
localized from itinerant polarons and to discern tightly bound polaron pairs (bipolarons) 
from simple polarons by such methods as photoemission spectroscopy. In section 5 we will 
study the time evolution of the charge redistribution and of the lattice deformation and 
discuss the renormalization of the intrinsic lattice vibrational frequency and of the electron 
hopping integral in the course of the charge transfer process. 



II. THE MODEL 

The smallest system on which polaronic features can be studied as far as self-trapping, 
localization-delocalization crossover, frequency renormalization and dephasing of the corre- 
lated charge-deformation dynamics is concerned, is the so called two-site polaron system. It 
consists of two elastically uncoupled adjacent diatomic molecules with atoms having mass M 
and an intrinsic bare frequency of intra- molecular vibrations ioq. Electrons can hop between 
those two molecules having a bare hopping rate t. In the Holstein model for such a system, 
the strength of the coupling constant of the density of charge carriers to the intra-molecular 
deformations is denoted by A. The model Hamiltonian for such a system is then given by 

a cr a 

+ Y(^i + ^l) + y^oK + ul)+U (nitnii + natna^) (1) 

where riiu = cl^c^^ denotes the density of charge carriers at molecular site i and having spin a. 
The intra-molecular deformations are denoted by ui and U2 respectively. We assume a simple 
form of repulsive intra-molecular Coulomb forces of strength U. The above Hamiltonian 
Eq.(|l]) has been extensively studied by a number of autors |3|-[T0| who, by means of exact 



diagonalization, have examined various aspects of the polaron problem. Diagonahzation 
of the above Hamiltonian can be rendered more efficient when decomposing it into a 
term containing the symmetric in-phase lattice vibrations characterized by a wavevector 
g = and antisymmetric out-of-phase vibrations characterized by a wavevector q = tt. The 
Hamiltonian then can be separated into two independent contributions H = Hx + Hy, given 
by: 



(2) 
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and a^'^^ and 6^^-* denote the annihilation (creation) operators of the quantized lattice fluctu- 
ations with momentum q = tt and q = respectively. Since in Hy the phonons couple only 
to the total charge n = ^ia of the system, Hy can be diagonalized exactly since it repre- 
sents simply a shifted oscillator, corresponding to the two diatomic molecules of our system 
having their equilibrium positions shifted by equal amounts = {Xn)/2MujQ). Hx on the 
contrary contains the full dynamics of the system and has to be diagonalized numerically in 
terms of a judiciously chosen set of basis states (see ref. such as 



\Ix,Iy) 



Nx,Ny 

oo ^ ^ 
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for the one-electron two-site problem. \Nx,y) denote the eigenstates of the oscillator part 
of Hx,Y for A = 0. The coefficients (3iy,Ny known analytically from the expansion of a 



shifted oscilator states in terms of the excited states of the unshifted one. For /y = we have 
I3o,Ny = ~ {cP' /2)a^"^ / \/Ny\- The coefficients a^^J^ are determined by diagonalyzing 
numerically Hx with in a truncated Hilbertspace of states with up to 100 phonon states, 
i.e., < Nx < 100. The procedure of separating off in the Hamiltonian a part which can 
be diagonalized exactly represents a substantial reduction in the numerical work of studying 
polarons not only for our two-site model but in general W% - 



Applying the standard LF approach to this problem amounts to using a representation 
in which the molecules exist in a set of oscillator states: those not containing any electron 
and being described by the full set of oscillator states labeled by and 
those containing an electron on one of the two sites and being described by shifted oscillator 
states ± Xo))m. and \^{Y — yo))„, those containing two electrons on one of the two 

sites and being described by |$(X ± 2Xo))m and \^{Y — 2yo))n and finally those containing 
one electron on each site and being described by |$(X))m and \^{Y — 2Fo))n- The shift in 
equilibrium positions for the two oscillators are given by Xq = Yq = Xn/ \/2MuJq and the 
shifted oscillator states are defined by 

\^X - Xo))™ = ^——J—exv{a\a^ - a])|$(X))o 
\Jm\ 



= ^^^4^E^exp(-aV2)|$(X)), (7) 

In such a representation the Hamiltonian Hx becomes 

Hx = tE(ni^ + n2a) - tJ2ic\^exp{a[a'' - a\)c2^ + H.c.) 

-Spn - {2ep -U)Y Ui^Uii + fiujQ{a)a + ^) (8) 

i 

in which the intricate dynamics of the polaron problem is contained in the modified hopping 
term (the second term in Eq.(^)) showing a concomitant transfer of charge and deformation. 
The standard LF approach to this problem is then to consider, to within a first approxi- 
mation, that the deformation of the molecules follows instantaneously the motion of the 
electrons, that is to say without any emission or absorbtion of phonons during the process 



7 



which transfers an electron from one site to the other. For our two-site polaron model this 
implies states of the form 

\LF)t. = ^ {clmX - Xo))m ± clmx + Xo))m) MY - Yo))n (9) 

which amounts to disentangle the correlated hopping term in the form 

Ci^exp{a[a^ — a])c2o- ^ 
cLc2<x(ea;p(a[a"^ - a\))n^^ = tlpc\„C2„ (10) 

This way one obtains an effective Hamiltonian to start with and subsequently treats the 
remainder of the full Hamiltonian in a perturbative way - corresponding to an expansion 
in terms of 1/A [|12|. The average in Eq.(|l^) is taken over the free phonon Hamiltonian 
Hph, given by the last term in Eq.(^ and which leads to t^F — texp{—2a'^). For a many 
polaron problem, polaronic states may become unstable with respect to bi-polaron forma- 
tion provided the Coulomb repulsion U is overcome. This then leads to states of the form 
^^^^^(X ± 2Xo))m|'^'(^ ± 2Fo))n- Upon eliminating single polaron states in the Hamilto- 
nian, Eq(|T|) one derives [|T^ an effective Hamiltonian describing bi-polaronic hopping with 
a hopping integral t^p = [t"^ / {2ep)exp{—Aa'^) for U = 0. 



III. INTERSITE CHARGE TRANSFER AND MULTIPLE HOPPING 

PROCESSES 

It has remained a question of dispute as how good the LF approach is and whether 
in the extreme polaronic limit the LF approximation, i.e. the substitution described by 



Eq.(lO), becomes exact. If this is not the case, then a perturbation theory in terms of 1/A 
can not be applied to the polaron problem. It has for a long time been taken as granted 
that the problem of a single polaron in an infinite lattice of finite dimension is described by 
polaron band states having a k dependent dispersion identical to that of the bare electron 
but with a reduced bandwidth 2zt*, where z denotes the coordination number. There are 
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now indications from exact diagonalization studies of finite clusters that this is not 
correct and that the k dependence of the dispersion differs significantely from that of the 
bare electrons. From exact diagonalisation studies P,p^ of the two-site polaron problem we 
know that over a large regime of parameters a and 7 the LF approximation gives reliable 
results as far as the energies of the low lying excitations are concerned. Nevertheless there 
can be serious discrepancies between the LF approach and the exact results when considering 
the eigenstates of the polaron problem. Let us illustrate theses discrepancies in such physical 
quantities as 

i) the kinetic energy of the electrons, 

ii) the occupation number of electrons as a function of wave vector, 

iii) the wave vector dependence of the spectral function measurable by angle resolved 
photoemission experiments, 

iv) the importance of renormalization of the phonon frequencies and of the electron 
hopping integral 

v) the retardation between the dynamics of the charge carriers and of the lattice defor- 
mation which accompanies them. 

Let us start by considering the static correlation function te///2t = —Ekin/2t = (c|^C2o.) 
which describes the kinetic energy in units of 2t and which in the extreme polaronic LF limit 
(Eq.([l0|)) becomes t\p/t for single polaron hopping and t*^p/2ept if it is bipolarons which 
hop. It is for this reason that sometimes this correlation function is associated with some 
effective hopping integral teff/2t which also in the limit a ^ is physically meaningfuU 
since it tends to the free electron value equal to \. Apart from these extreme limits, we shall 
see below, the interpretation of this correlation function in terms of an effective hopping 
integral has to be modified. Evaluating this correlation function for the two-site polaron 
problem by exact diagonalization of the Hamiltonian Eq.(|l|) we notice that, as a function of 
a and 7, the exact result for Ekin/2t in general strongly deviates from the extreme polaronic 
LF limiting behavior. We can hence not expect that the LF approximation is applicable 
even in conjunction with perturbative corrections in terms of 1/A. 



In Fig. (la) we plot teff/2t as a function of a for different adiabaticity parameters 7. At a 
first glance we get the impression that upon going towards the extreme anti-adiabatic limit 
i.e. 7 — > and for a fixed value of a we approach the extreme polaronic LF behavior. Yet as 
can bee seen from Fig. (lb) and contrary to a well established belief on this matter, the LF 
result is approached slower and slower for a fixed 7 ^ 1 as a gets bigger and bigger. These 
discrepancies with the LF results are not only quantitative but are qualitative in nature as 
shows the 7 dependence of te///2t. The reason for this qualitatively different behavior lies 
in the LF approach itself which considers the electrons to be in localized states of the form 
given by Eq.(^ in which practically the entire charge of the electron and entire deformation 
of the molecules remains restricted to the molecular site on which the electron sits. This 
leads to tefj/t = t}^p/t = q{^{X — Xq)\^{X + Xo))q at zero temperature which noticeably 
differs from 

t,jf/t= +(vl/(X)|M/(X))o (11) 

where |\1/(X))^ denote the exact eigenstates which replace \^{X ±X^))m in the expressions 
for the LF approximated eigenstates given by Eq.(|^). As can be seen from Fig. (2) the exact 
eigenfunctions \1/^(X), corresponding to the exact eigenstates |\E'(X))^, differs significantly 

from the LF approximated eigenfunctions ^^(X) which for the groundstate becomes 

1 

^oiX) = {^^y exp[-{X ±XofiMu;o/n)] (12) 

In contrast to $^(X) the exact eigenfunctions \l/o (X) show a substantial deformation of 
the molecule adjacent to the one where the electron actually sits. This is born out in the 
smaller of the two peaks of \l/o (X) for negative values of ^ = X ^ Muo/h ~ —Xq^ Muq/H. 
It is the presence of this second peak which gives rise to a substantial increase in the value 
of teff/t over the LF approximated t}^p/t. As the coupling strength increases the value of 
Xq increases roughly linearly with a. This leads to a separation of the two main peaks of 
\E',|(X) and (X) respectively, which results in an exponentially small overlap of them. 
On the contrary the overlap of the main peak of \E'(J'(X) with the secondary peak of \E'o (X) 
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is of order unity since the weight of the secondary peak depends only weakly on the value of 
a. Since this second contribution to ^e///^ always remains much bigger than the first one, 
we can thus never hope to recuperate the LF result t\pjt for whatever fixed value of 7 < 1. 

The behavior of the oscillator wavefunction expresses a dynamical delocalization of the 
polaronic state which has been realized a long time ago and for which variational calculations 
have given rather satisfactory results for the groundstate (see for instance ref From the 
above study of te///2t and of the oscillator wave function it becomes clear that in 

general effects of dynamical delocalization of the electron cannot be obtained by perturbative 
expansions in terms of 1/A around the LF approximated oscillator wave function, even in 
the extreme antiadiabatic limit 7^1. Our findings suggest that it is energetically favorable 
for the polaron to be partially delocalized and to transfer its charge by multiple hopping 
processes because it is that which gives rise to the monotonically increasing behavior of 
teff It as a function of 7. These multiple hopping processes play, as we shall see in section 5, 
an important role in the correlated charge - deformation dynamics of the polaron problem. 

As concerns tejf/t for the two electron two-site system we notice from Fig. (3a) a behavior 
similar to that for the one electron problem. However as can be seen from Fig. (3b), in 
contrast to the one electron problem, for two electrons it varies linearly with 7 which is an 
indication that it is bipolarons rather than polarons which hop. As 7 increases, teff/t tends 
towards a constant equal to | which is indicative of uncorrelated hopping of the two electrons 
in the system. There is a significant difference in slope of this correlation function calculated 
exactly and of that determined by the LF approach. Again, for reasonable parameters 
where small bipolarons are stable (such as a = 1.2 and 7 ~ 1) we find that the exactly 
calculated correlaton function is orders of magnitude bigger than its LF expression for the 
same parameters. 

In concluding this section we want to point out that although polarons as well as bi- 
polarons exist over a large regime of parameters a and 7, their description in terms of 
the standard LF approach is generally invalid not only in the physically interesting regime 
(which lies outside the extreme strong coupling antiadiabatic limit 7^1) but particularly 
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in that limit, for which it largely overestimates the degree of self-trapping. The failure of 
calculating tf.ff/2t qualitatively correctly within the LF approach was already recognized 
by Fehske et al. |T5[ who attributed it to the zero-phonon approximation, inherent in the 



substitution Eq.(lO). They showed [0] that the distribution of the weights of the N phonon 
state in the ground state shift its maximum to higher values of with increasing values 
of Ep/t = 2o?^. For a fixed values of a the maximum of this distribution shifts to higher 
and higher value of as 7 decreases, thus rendering the zero-phonon approximation less 
and less justified. In order to visualize this behavior we plot this distribution Pin) = 
J2Nx=oi\'^oNx\'^ + I'^cwvx Hl-^o^-^x P a fuuctiou of for a set of values of a (Fig. 4(a)) 
and of 7 (fig. 4(b)). From Fig.4(b) we can see from a comparison of the LF approximated 
distribution with the exact one for a particular choice of a = 1.2 how the LF aproach 
overestimates the the self-trapping, with a maximum in P{N) occuring at a significantly 
higher value of than is the case in the exactly calculated P{N). 

The findings of this section moreover suggest that the effect of dynamical delocalization 
of the charge carriers in polaronic systems leads to a strengthening the phase coherence 
of polaronic states. If this is so, then we expect that for a given set of parameters a 
and 7, states described by the LF approach become unstable, i.e., have their spatial phase 
coherence destroyed as the temperature increases, while the exact states will maintain this 
phase coherence up to much higher temperatures. The key to this question lies in the wave 
vector dependence of the one particle spectral function which will be discussed in the next 
section and of the occupation number n„{k) = {c\^c^J) of the electronic charge carriers. 
Evaluating ncr{k) for the two lowest eigenstates within the LF approach, that is to say 
with respect to the two eigenstates given in Eq.(|^) with m = n = we obtain na{k = 
0)+ = n„{k = 7r)_ = i(l + exp(— a^)) and n„{k = 7r)_ = n^{k = 0)+ = |(1 — exp{—a'^)) 
respectively. This inversion of the occupation numbers between the ground state and the 
first exited state leads to the result that at low temperatures (bigger than the difference in 
energy of the two lowest levels but small compared to the phonon frequeny uq) n{k) = i, 
independent on the wavevector k and thus identical to the result for localized polarons, i.e. 
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for t = 0. Evaluating ncr{k = 0, 7r)-|- exactly we see from Fig. (5 ) that the exact result differs 
from that of the LF approach qualitatively with n„{k = 0)± not only being larger than 
rif^i^k = 7i)± for the ground state but also for the first excited state, provided a > 0^(7) 
with acril) increasing monotonically as 7 decreases. 0^(7) is determined by that value of a 
for which teff/t in Fig. 1(b) tends to the saturation value |, given by the free electron limit, 
for a given value of 7. It is this fact which gives the polaron a stronger dynamical coherence 
than what we would expect on the basis of the LF approach and which shows up in the 
temperature variation of the effective polaron hopping integral te//, illustrated in Fig. (6). 
We notice that in the weak coupling regime (a <^ 1), where the electrons behave as quasi-free 
charge carriers, tg// decreases with increasing temperature. This is precisely what we expect 
if the phase coherence of the electron is destroyed by thermal fluctuations. On the contrary, 
in the polaronic regime (a > 1) we notice an initial small increase of tg// with increasing 
temperature (corresponding to a decrease in the kinetic energy of the electrons!) which 
suggests that the dynamical coherence of the polaron increases with increasing temperature. 
This holds true for temperatures up to some characteristic temperature above which this 
coherence is definitely destroyed resulting in a decrease of t^jf with increasing temperature. 
This once more strongly contrasts with the result obtained within a LF approach which 
treats the polarons in terms of band states and which consequently leads to the well known 
reduction of t*^p = exp{— 20^ coth{(3ujQ/ 2)) with increasing temperature (see Fig. (6)). 

These features are of course quite general and will remain when, instead of studying the 
two-site polaron problem, one deals with the polaron problem in an infinite lattice of finite 
dimension. There the LF approach is in fact known to give rise to results which contradict 
our findings for the two-site polaron problem such as the dependence of the occupation 
number on the wavevector which, apart from a small step at the Fermi vector, is a fiat 



function throughout the Brillouin zone for zero temperature [0. For any realistic finite 
temperature, which clearly would be larger than the difference in energy of the two lowest 
lying eigenstates of the Hamiltonian Eq.([I|) (i.e ~ 2texp{—2a^)) , the LF approach thus would 
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lead one to expect n„{k) to be equal to |, independent on k and not to show the slightest 
anomaly at the Fermi wavevector. Furthermore as concerns the temperature variation of 
teff discussed above for the two-site polaron problem one expects also for an infinite lattice 
a mobility which increases with increasing temperature. The opposite behavior is found 
in the classical works on that issue and being based on the LF 1/A perturbative approach 
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An important feature in polaron physics is the crossover from small polarons in the anti- 
adiabatic strong coupling limit to large polarons in the adiabatic strong coupling limit. This 
crossover has been studied in great detail and can be seen most clearly if t^ff is plotted as a 
function of 6p/t = 2a'^ujQ/t [jl^ rather than as here as a function of a. It turns out that the 
crossover is rather abrupt in the adiabatic strong coupling limit and becomes more spread 
out when going to the anti-adiabatic limit. 

In the following two sections we shall discuss finer details of polaron dynamics which 
should show up in the spectral properties and the correlated dynamics of the charge carriers 
and of the molecular deformations which accompany them. 



IV. ITINERANT VERSUS LOCALIZED POLARONS 

One of the key questions in the physics of small polarons is how to discern itinerant 
from localized polarons. In principle the answer to this question should be contained in 
the one particle density of states, which can be measured by photoemission experiments. 
Considering this problem within the LF approximation, the polarons are described by band 
states which are characterized by the total momentum k = k + q, where k denotes the 
momentum of the electron and q that of the deformation which accompanies it. As an 
illustrative example for the one electron two-site polaron problem this implies that the LF 
groundstate \LF)qq = |LF)k=o, given by Eq.(D, can be written as 

\LF).=o = ^ (4|$(x))o + 4 - ro))o (i3) 

where 
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c[o,,), = (cL ± cL) , |$(X))o,. = ^ mX - Xo))o ± |$(X + Xo))o) (14) 
Hence the scattering cross-sections 

1 r 

^Sval,<x2,...a^(^.^) = J dT exp(i^r) ^(cL(r)Cfc^(0))7v^(r) 

^S^:l,.„....,(fc,^) = ^"^^ / dr exp (z^r) ^(c,,(r)cL(0)),v^(r) (15) 

measure the intensity of emitted and absorbed electrons with momentum k and spin cr in a 
polaronic state containing electrons with spins (Xi, cr2, ...ctat in the groundstate. 

Let us first consider the case of an isolated polaronic center with one and respectively 
two electrons present. We shall from now on consider only the case U = 0. The two 
cross-sections, given in Eq.(0), are then exactly determined by 

„2« 



^ ',1T ^ 2^ex77|^-« oi^u; - fcp -h lujq) 

„2; 



^n;2n(^) =Eea:p(-«')^5(^-3£,T/^o) (16) 

which gives the well known Poisson distribution of the phonon modes locked together in the 
construction of the coherent Glauber states which define localized polarons. We shall now 
show that the overall features of these spectral functions are preserved when considering 
a system of itinerant polarons. For that purpose let us consider the spectral function for 
electron emission from a one polaron ground state in the two-site polaron problem within 
the LF approach. We then obtain 

I even I' 

oo oo ^2(i+/') 

A%\k = 7:,uj),= Y.Y.^xp{-2a^)——5{uj~e,-{l + l'U) (17) 

I odd V • 

if the system is in the ground state and y4''^|^(fc = 0,ti;)i (= y4*-^|''(A; = 11,00)0) and y4^^|^(fc = 
7r,ci;)i (= A\^ [k = 0,ti;)o) if it is in the first excited state. This shows that the bulk of 
the spectrum of itinerant polarons (i.e. t 7^ 0) is identical to that of localized ones except 
for the lowest energy part of the spectrum i.e. for small values of /. For that part of the 
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spectrum the spectral weight ior k — in the ground state equals ~ exp — {2a^) while it is 
identicaly zero for k = n. The opposite is true for the first excited state. These are precisely 
the features expected for quasi-particles whose weights for their coherent part is strongly 
reduced. The exact results of the spectral functions bare this out and in Fig. (7a) we present 
, ct;)o for values of 7 and a which characterize well defined polaronic states. In 
order to compare A^f^\k = 0, a;) with ^'•^■'(a;) for localized polarons (i.e. t = 0) as well as 
with {k — TT, u) we do this for finite rather than zero temperature since for all intents 
and purposes the level splitting between the fist two lowest energy levels (corresponding to 
the bandwidth 2t* of the coherent polaron motion) is extremely small, since under normal 
conditions we are dealing with temperatures T such that 2t* <^ ksT <g loq. Under those 
circumstances the LF approach yields scattering cross-sections A^f^\k,uj) which become 
wave vector independent and thus undistinguishable from purely localized polarons. The 
exactly calculated spectral functions become, up to a A; dependent scaling factor, practically 
identical with A{u!) for the entire frequency regime except for the low frequency part, where 
however the spectral weight is extremely small (see Fig. 7b). These findings suggest that 
small polarons might never exist in form of coherent Bloch like states, a feature which is 
supported by the dynamics of the polaron motion, which will be discussed in section V. 

One of the prime problems in the physics of the many polaron problem is the question 
of how to discern a system of essentially non-interacting polarons from one in which they 
sense a strong attraction between them and which ultimately leads to the formation of 
tightly bound pairs (bipolarons) . In principle this question can be resolved by photoelectron 
spectroscopy which measures the one particle density of states for polarons. In fact what is 
required in order to discriminate between a many body ground state containing essentially 
unpaired from one containing essentially paired polarons is both, photoemission and inverse 
photoemission spectroscopy. Testing the system with solely photoemisson we can not decide 
about the nature of the many body ground state. This can be illustrated on the basis of the 
two-site polaron system containing one and respectively two electrons. The scattering cross 
section looks essentially similar for the two cases. In Fig.(8,a-d) we plot for different values 
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of coupling strength a the photoemission spectrum A (A; = 0,0;) where an electron is 
emitted out of a two-electron two-site polaron system. This spectrum has the same form as 
,uj) (see Fig. (7a)) which corresponds to the situation of an electron emitted out 
of a one-electron two-site system. On the contrary, the inverse photoemission spectrum can 
clearly distinguish if the final state corresponds to a state with essentially uncoupled polarons 
or polaron pairs. Looking at this spectral function for the one-electron two-site system i.e., 
y4*-^|''(fc = 0,Ci;) we see from Fig.(9a-d) that, as we increase the coupling strength, a two 
peak structure emerges. The energy gap which separates this two peak structure amounts 
to 2ep, which represents precisely the binding energy between two polarons in the strong 
coupling limit. The low energy peak of this spectrum corresponds to a final state given by 
a biporaron in its bonding singlet state, having energy ~ ASp. The high energy peak arises 
from a bipolaronic state in its antibonding singlet and respectively triplet state having an 
energy ^ 2ep |jl8|. The difference in energy between these two contributions in the spectral 
function for inverse photoemission is hence just the binding energy of a bipolaron and thus 
can serve as a signature for a bipolaronic many polaron ground state. 

It is even moreover illustrative to study the one particle density of states for those 
situations. In Figs.(10,a-d) we plot the density of states per spin for the one electron two- 
site system at zero temperature 

pi(cu)= {A%\k,u), + A%\k,u),) (18) 

fc=0,7r 

as a function of coupling strength a and for a given fixed value of adiabaticity parameter 
7 = 1.1. For small values of a the density of states is characterized by two peaks centered 
around energies c<j ~ and uo ~ 2t^ which in essence represents quasi free electrons in the 
bonding and respectively anti-bonding states for this small system. As a increases these two 
peaks spread and eventually evolve into two well separated peaks. This density of states 
is similar to that of the Hubbard model in the dilute limit (i.e., far away from half filling) 
which shows indications for the presence of an upper and a lower Hubbard band, separated 
by the Hubbard U repulsion energy, but without any clear gap between those two bands. 
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Similarly to the Hubbard problem, the two peaks in the polaron problem are separated by 
the energy of attraction between two polarons, which in case of strong coupling (such as 
illustrated in Fig.(9d) for a = 2.2) amounts to 2Sp. 

A similar density of states per spin is obtained for the two electron two-site system 



which we illustrate in Figs.ll(a-d) as a function of coupling strength a and for a given fixed 
value of adiabaticity parameter 7 = 1.1. We again see a similar evolution of the density of 
states as we increase a but now the separation in energy is equal to iSp, which can be clearly 
distinguished from the case representative for uncoupled polarons (the one electron two site 
problem considered above) . A clear energy gap given by 2£p now appears separating the low 
and high frequency peak structures. Integrating p2{i^) up to some value uj = Sp such that it 
becomes equal to 2 (which corresponds to the half filled band case and to two electrons in 
our system) we find that lies precisely in the middle of the gap of the density of states. 
This again is reminiscent of the problem of strongly correlated electrons for the half filled 
band Hubbard model. These features examined here for the two-site polaron system are 
expected to hold true generally for any interacting many polaron system on a lattice. 



The polaron problem presents a highly non-linear dynamical system in which the charge 
and deformation fiuctuations are intricately coupled together. This leads to a dynamics 
of the molecular deformations being driven by the dynamics of the charge carriers in the 
strong coupling anti-adiabatic regime (a ^ 1,7 <S 1) . On the contrary it is the dynamics 
of the molecular deformations which drives the dynamics of the charge carriers in the weak 
coupling adiabatic regime (a <^ 1, 7 ^ 1). In general this leads to a dynamics for the charge 
and of the molecular deformations which is composed of a common slow oscillation and fast 
ones superposed on it. The fast oscillations for the charge dynamics have a frequency t of 




(19) 



fe=0,7r 



V. CORRELATED CHARGE-DEFORMATION DYNAMICS 
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the order of t, while those for the deformation dynamics have a frequency uq which is of the 
order of cuq. 

In this section we shall study the evolution of those dynamical properties of small po- 
larons when we go from the strong coupling anti-adiabatic regime to the strong coupling 
adiabatic one, with a special emphasis on the behavior in the crossover regime. In order to 
illustrate this behavior we evaluate the time dependent correlation functions for the charge 
redistribution and molecular deformations: 

Xnn(r) = (Mr) - n2a(r))M0) - n2.(0))), xM = {X(t)X{0)) (20) 

In Figs.l2(a-d) we plot these correlation functions (normahzed with respect to their r — 
values) as a function of time r in units of cuq for different values of the 7 and fixed a. We 
notice that the charge dynamics qualitatively tracks globally the behavior expected on the 
basis of the LF approximation in the anti-adiabatic hmit (i.e., for 7 = 0.1, Fig.l2(a)) but 
with superposed small amplitude fast charge oscillations wih a frequency t which is large 
compared to the unrenormalized electron hopping integral t. The dynamics of the molecular 
deformation follows in a coherent fashion that of the charge and exhibits superposed molec- 
ular vibrations with a frequency cuq ~ cuq. As we increase 7 (i.e., for 7 = 1.1, Fig.l2(b)), the 
amplitude of the fast charge oscillations become increasingly important and their oscillation 
frequency t decreases. The opposite behavior is obtained for the molecular deformation 
oscillations, whose frequency ujq of the fast oscillatory behavior increases while the corre- 
sponding amplitude diminishes. The system is then no longer described by the LF approach 
(as we can see from the comparison made for this case) and the charge dynamics is now 
controlled by multiple hopping processes with concomitantly reduced amplitude fluctuations 
of the molecular deformations. Upon further increasing 7 we arrive at a situation where the 
frequencies and amplitudes of those two dynamical variables, characterizing the charge and 
deformation, become comparable to each other i.e., for 2t ~ ujq and we enter the crossover 
regime (i.e., for 7 = 1.6, Fig. 12(c)). This regime is characterized by a temporarily alternat- 
ing behavior between essentially self-trapped anti-adiabatic small polarons (manifest in a 
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substantial reduction of the amplitude fluctuations of the charge dynamics) and a behavior 
reminiscent of itinerant adiabatic polarons as seen upon a further increase of 7 (ie.,7 = 2.0, 
Fig. 12(d)). Such fluctuations in the amplitudes of the charge occur over a time scale which 
is large compared to the inverse hopping rate t and is accompanied by phase slips in the fast 
oscillatory behavior of the charge as well as the deformation fluctuations. 

The evolution of the correlation functions for the charge and deformation fluctuations of 
the two-electron two-site problem as a function of the adiabaticity parameter 7 and flxed 
coupling constant a follows a similar behavior to that of the one-electron two-site problem, 
as shown in Figs.l3(a-d). Again we can identify a crossover between small bi-polarons and 
essentially uncorrelated two electrons which occurs for 7 ^ 1.3 for the particular choice 
a = 0.6 

This crossover between self-trapped polarons, respectively bi-polarons, and quasifree 
electrons in the phase space of a and 7 corresponds to the characteristic value of 7 where the 
kinetic energy of the electrons E^in = —2t{c\ 102,1) approaches its maximal value of the free 
electron limit, while the potential energy of the electrons Epot = — 2A((?2i,| + ni Jmi) tends 
to its minimum value — ^Ayo((^i,T + ^^1,1)); obtained in the limit 7 ^ 00. This can be seen 
from the behavior of Ekin illustrated in Fig. 1(b) and of Epot depicted in Fig. 14 for a = 1.2 
and 7 ~ 1.6 for the one-electron two-site problem. The crossover between self-trapped small 
polarons and quasi-free electrons thus appears to be driven by a competition between the 
kinetic and potential energy of the electrons; the flrst one favoring a delocalization of them 
while the second one inciting them to localize on the molecular sites. In an inflnite solid 
state system such a scenario would suggest a quantum phase transition between a metal and 
a polaronic insulator as proposed a long time ago by Landau and Froehlich ||20|. For systems 



with low carrier concentrations such localized polarons have been verifled experimentally in 



metal halide where optically excited excitons get localized |2T[]. There is at present no exact 
theorem as to whether pure electron phonon-systems can show such a polaronic insulator. 
The present exact proofs against that [^] hinge on suppositions of the phonon spectra which 
may not be realistic for real materials. 
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We finally should like to point point out that no significant changes in the phonon 
distribution of the displaced oscillators are observed when going through this crossover 
regime. This can be verified from the plot in Fig.4(b) of P{N) for a fixed a = 1.2 and upon 
varying 7. 

In Table I we summarize our findings on the charge-deformation dynamics of a pola- 
ronic system. We compare for that purpose the renormalized vibrational frequency ujq, the 
renormalized electron hopping integral t, the kinetic energy of the electrons teff/t and the 
physical charge transfer rate t* for different values of the adiabaticity parameter 7 and for 
a fixed coupling constant a. In the limit of strong adiabaticity, t* tends to the LF value t^F 
and follows as a function of 7 the behavior of {Ei — Eq)/2 which denotes the difference in 
energy of the two lowest eigenstates. We also indicate the spectral weight Z of the lowest en- 
ergy contribution to the scattering cross-section and notice that for the anti-adiabatic limit 
it scales with t* /t. This is an indication that the low frequency part of the scattering cross- 
section corresponds to coherent states but with a spectral weight which is extremely small. 
It is presently not clear whether such weak coherent features, characteristic of itinerant small 
polarons, will persist if one treats the polaron problem on a infinite lattice. Calculations 
based on infinite dimensions EM which show such itinerant behavior do neglect totally the 



frequency renormalization of the vibrational motion of the atoms, which we consider as the 
prime cause for dephasing of the correlated charge-deformation dynamics and ultimately 
believed to be responsible for the destruction of itinerant polaronic states. 

The strong coupling between the charge and the molecular deformations thus manifests 
itself not only in a strong renormalization of the molecular vibrational frequency Uq becoming 
Uq but also in a strong renormalization of the intrinsic hopping integral t renormalized into 
t. As we go from the anti-adiabatic limit toward the adiabatic one (for fixed value of 
a), we observe a substantial decrease in t/t and a concomitant increase in ujq. These are 
effects which should be observable by spectroscopic measurements such as infrared or Raman 
scattering for the vibrational modes. 

^From inspection of Figs. (12, 13) we notice the sizeable dynamical delocalization of the 
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polaron and bi-polaron respectively as we approach the crossover regime. In the crossover 
regime itself this delocalization alternates between partly quasi-static delocalization, suggest- 
ing almost localized yet extended polaron states and dynamical delocalization reminiscent 
of almost free carriers, which however remain dynamically tied to a given molecule for some 
appreciable time. To be more specific, for the particular case illustrated in Fig. 12(c), we 
find for the period of time when a quasi-static polaron is stable, a charge distribution given 
by: (^i,(t) — 0.75 and (^2,0-) — 0.25. On the other hand for the dynamically delocalized 
polaron we find that the charge distribution fluctuates over a characteristic time given by t 
between {ni^„) ~ 1.0, {n2^a) — and(ni^o-) — 0.5, (?t.2,o-) — 0.5. Such temporal fluctuations 



were initially hypothezised by us a long time ago [0 which led to the Boson- Fermion model 
for intermediary coupling electron phonon systems, which may have some relevance for our 
understanding of High Tc superconductors . 



VI. SUMMARY 

The main objective of this work was the study of the intricate dynamics of the polaron 
problem involving the dynamical behavior of the charge carriers and that of the molecular 
deformations which surround them. We find that in the anti-adiabatic regime for small 
polarons the molecular deformations follow in a coherent fashion the redistribution of the 
charge, while in the adiabatic regime it is the charge redistribution which follows the molec- 
ular deformations. In the crossover regime between those two limiting cases we find that 
the dynamical behavior of the polaronic charge carriers alternates between self-trapped po- 
larons and almost free carrier behavior. The time scale over which these different behaviors 
are realized is typically an order of magnitude bigger than the intrinsic hopping rate i.e., 
of the order 10 x 27i/t. This crossover regime is characterized by strong renormalization of 
the intrinsic hopping integral as well as of the bare phonon frequency, which in this regime 
become equal. Phase slips in the fast oscillatory components of the charge and molecular 
deformation fluctuations are the result of this. Such effects are expected to be essential 
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for a proper description of polaron damping, sofar having been treated only within the LF 
approach and which is unable to account for the effects described here. 

The question of how to distinguish itinerant polarons from localized ones was studied 
here on the basis of the one particle spectral function and its temperature dependence. We 
showed that at zero temperature the respective spectral functions for localized and itinerant 
polarons differ from each other only very slightly except for the low frequency regime of the 
spectral function, where they show an increased spectral weight for small wave vectors if 
the polarons are itinerant. As the temperature increases this difference disappears and it 
is, in principle, no longer possible by spectroscopic means to distiguish between localized 
and itinerant small polarons. This may explain the puzzling results in the photoemission 
spectra for certain High superconducting cuprates for which a wave vector independent 
spectral function was observed in the normal state and consequently was interpretated as 



indications for localized charge carriers [0. The question of a polaronic insulator versus 
a polaronic metal has been touched upon here only from the point of view of the single 
particle properties. The polaron problem is however presents a problem of electrons in a 
system with impurity centers with dynamically varying energies and thus contains features 
similar to those of the Anderson localization. The relevant quantity to be studied hence 
is the conductivity. So far a few attempts in this direction have been made on the basis of 
exact diagonalization studies in finite systems attempting to determine whether there 
is or not a finite Drude component in the optical conductivity. 

Finally, our exact diagonalization studies on the two-site molecular Holstein polaron 
model permitted us to discuss the limitations of the standard LF approach. Our rather un- 
expected and perhaps widely unrecognized findings are that this approach, which is generally 
believed to become exact in the limit of anti-adiabaticity and an electron phonon coupling 
going to infinity, actually diverges most from the exact results precisely in this limit. The 
reasons for that can be traced back to the zero-phonon approximation inherent in the LF 
approach, based on the relation Eq. (10) and which, with increasing coupling strength, is 
increasingly strongly violated flSl . 
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Our analysis of the various spectral functions and the density of states shows that the 
major part of the spectrum must be considered as being due to incoherent rather than 
coherent polaron dynamics; the latter having vanishingly small spectral weight of order 
exp(— 2q;^). This result confirms our earlier findings on the many polaron problem for 
infinite lattices and examinations of the single polaron problem in infinite dimensions 
2|. 
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FIGURES 



Fig.l te///2t = —Ekin/2t = (c|^C2o-) for a single electron as a function of a for different 
adiabaticity parameters 7 for the two-site polaron problem (a). Comparison of tefj/t eval- 
uated by exact diagonalization and its approximative value given by the LP approach as a 
function of 7 for several values of a (b). 



Fig. 2 The exact oscillator wave function ^,[(0 with ^ = X^Muq/Ti for a single electron 
in the two-site polaron system for different values of a. 

Fig. 3 teff/t for two electrons as a function of a for different adiabaticity parameters 7 for 
the two-site polaron problem (a). Comparison of teff/t evaluated by exact diagonalization 
with its approximated value given by the LF approach as a function of 7 for different values 
of a (b). 

Fig. 4 (a) The phonon number distribution P{N) for a given value of a = 1.2 and for 
various values of 7. Notice that as 7 decreases, the distribution function tends to that 
of localized polarons, practically indistinguishable from that of 7 = 0.1. (b) The phonon 
number distribution for a given value of 7 = 1.6 and for various values of a. Notice the 
in general significant difference between the LF approximated distribution function and the 
exact results, illustrated here for a = 1.2. 

Fig. 5 The wave vector dependence of the polaron occupation number n^°\k)± for the 
two lowest lying states. For comparison we also plot the result for the LF approximated 
states given in Eq.(|^). 
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Fig. 6 teff/t (normalized to its zero temperature value) as a function of temperature 
(in units of uiq) for fixed 7 = 1.1 and for different values of a and compared with the LF 
approach results. 

Fig. 7 The scattering cross-section {k = 0, a;)o for 7 = 1.1 and a — 2.2 and for ksT — 
O.la;o (a). Comparison of the scattering cross-sections A^f^\k = 0, a;)o and A^f^\7r,u!)o with 
each other as well as with that for localized polarons i.e., 74^^^(a;)o for ksT — O.luo and 
T = (b). These scattering cross-sections have been obtained by broadening the set of S 
functions by Gaussians of width Aa; = 0.1a;o. 

Fig.8(a-d) The scattering cross-section 74^2||(^ = 0,a;)o for photoemission from a two 
electron two-site system as a function of increasing electron phonon coupling a. These scat- 
tering cross-sections have been obtained by broadening the set of 5 functions by Gaussians 
of width Acu — OAcuq. 

Fig. 9 The scattering cross-section 74'-^^ (A; = 0, a;)o for inverse photoemission from a 
one-electron two-site system as a function of increasing electron phonon coupling a. The 
scattering cross-section has been obtained by broadening the set of 6 functions by Gaussians 
of width Acu — 0.1a;o. 

Fig.lO(a-d) The evolution of the density of states Pi{lv) for the one-electron two-site 
polaron problem as the electron phonon coupling a increases and for a fixed adiabaticity 
parameter 7 = 1.1. The densities of states have been obtained by broadening the set of 5 
functions by Gaussians of width Aa; = 0.1a;o. 
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Fig.ll(a-d) The evolution of the density of states p2{oj) for the two-electron two-site 
polaron problem as the electron phonon coupling a increases and for a fixed adiabaticity 
parameter 7 = 1.1. The densities of states have been obtained by broadening the set of 8 
functions by Gaussians of width Ao; = 0.1a;o- 

Fig. 12a Xnn and XxxifyK^"^) for a = 1.2 and 7 = 0.1 for the one-electron two-site 
polaron system. 

Fig. 12b Xnn and Xxx{t)/{^^) for a = 1.2 and tj = 1.1 for the one-electron two-site 
polaron system. 

Fig. 12c Xnn and Xxx{t)/{^'^) for « = 1-2 and 7 = 1.6 for the one-electron two-site 
polaron system. 

Fig.l2d Xnn and Xxx{'t)/{^^) for <^ = 1-2 and 7 = 2.0 for the one-electron two-site 
polaron system. 

Fig. 13a Xnn and Xxx{t) / i^"^) for a = 0.6 and 7 = 0.2 for the two-electron two-site 
polaron system. 

Fig. 13b Xnn and Xxx{'t)/{^^) for = 0.6 and 7 = 0.8 for the two-electron two-site 
polaron system. 

Fig. 13c Xnn and 

Xxxif)/ (^^) for OL — 0.6 and 7 — 1.3 for the two-electron two-site 

polaron system. 
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Fig. 13d Xnn and Xxx{t)/{X'^) for a — 0.6 and 7 = 1.6 for the two-electron two-site 
polaron system. 

Fig.l4 Comparison of the electronic kinetic and potential energy Epot — — 2A((ni,-f -|- 
ni^l)ui) ( as a function of 7 for a fixed a — 1.2 for the one-electron two-site problem. The 
crossover regime between self-trapped polarons and quasi-free electrons occurs when the 
kinetic energy tends to its maximal value while that of the potential energy tends to its 
minimal value. 
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TABLES 



1 = l/^u 


/// 








AE/21 


/*// 


Z 


0.1 


- 


1.03 


- 


1.00 10-1 


5.65 10-2 


5.61 10-2 


6.11 10-2 


0.3 


5.66 


1.08 


1.57 


1.90 10-1 


5.73 10-2 


5.66 10-2 


7.25 10-2 


0.5 


3.05 


1.14 


1.34 


2.82 10-1 


5.95 10-2 


5.98 10-2 


8.59 10-2 


1.1 


1.56 


1.42 


1.21 


5.64 10-1 


6.98 10-2 


6.79 10-2 


1.34 10-1 


1.3 


1.41 


1.52 


1.07 


6.20 10-1 


7.63 10-2 


6.90 10-2 


1.49 10-1 


1.6 


1.25 






7.12 10-1 


7.84 10-2 


7.85 10-2 


1.70 10-1 


1.7 


1.28 






7.38 10-1 


7.97 10-2 


8.03 10-2 


1.75 10-1 


2.0 


1.19 






8.00 10-1 


8.26 10-2 


8.26 10-2 


1.90 10-1 



TABLE I. The variation of the renormaUzed frequency of the deformation oscillations lvq and 
of the renomalized hopping rate t as a function of the adiabaticity parameter 7 for fixed coupling 
constant a = 1.2. Notice that as we approach the crossover regime, the time scales of these two 
oscillations become equal i.e., i/Coo 1. We also compare the frequency t* of the slow polaronic 
motion with the splitting of the two lowest eigenvalues AE/2t 

and teff (the electron kinetic energy). Notice that the spectral weight Z of the lowest 
frequency pole of the electron Green's function scales fairly well with the renormalization 
factor for the polaron bandwidth t*/t. 
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